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Abstract 

A method is proposed for producing monoergetic, high-quality ion beams in vacuum, via direct acceleration 
by the electromagnetic field of two counterpropagating, variable-frequency lasers: ions are trapped and accelerated 
by a beat-wave structure with variable phase velocity, allowing for fine control over the energy and the charge 
of the beam via tuning of the frequency variation. The physical mechanism is described with a one-dimensional 
theory, providing the general conditions for trapping and scaling laws for the relevant features of the ion beam. 
Two-dimensional, electromagnetic particle-in-cell simulations, in which hydrogen gas is considered as an ion source, 
confirm the validity and the robustness of the method. 

I. Introduction 

The recent developments in the technology of ultraintense infrared (IR) lasers have opened the way to compact 
particle accelerators for ions and electrons. In the case of electrons, the possibility of achieving relativistic radiation 
intensities (/ > 10 18 W/cm 2 ) motivated a number of proposals for direct acceleration in vacuum by the electro- 
magnetic (EM) field of one or more laser pulses [1] — [13]. However, indirect processes, in which the acceleration is 
obtained via laser-driven plasma waves [14], were shown to be more effective and reliable [15]— [18]. In the case of 
ions, for which relativistic radiation intensities (/ > 10 24 W/cm 2 ) are still beyond the limits of current technology, 
direct acceleration has been mostly unexplored (only very recently, direct acceleration by tightly focused, radially 
polarized lasers has been investigated [19]). Most of the proposals for laser-based ion acceleration rely on indirect 
processes, in which the ions are accelerated by the space-charge field in laser-irradiated solid targets (cf. [20]-[29] 
and, for a review, [30]-[32]), resorting to various physical mechanisms, such as, with increasing laser intensity, 
plasma-expansions [33]-[35], electrostatic shocks [36], [37], and the laser-piston regime [38]. In such processes, 
the properties of the accelerated ions are often difficult to control, making the production of high-quality beams a 
challenging task. In the spirit of conventional accelerators, better results would be obtained using schemes in which 
the acceleration is provided by a controllable wave structure. To accelerate ions in the nonrelativistic regime, a slow 
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wave is required, whose velocity increases as the ions accelerate. Possible methods have been proposed to drive such 
waves in plasmas [39]-[41]. Standing or slow waves can also be generated in the beating of counterpropagating 
lasers with equal or slightly different frequencies: a typical example is the optical injection of electrons into a 
plasma wave in a laser-wakefield accelerator (LWFA) [18], [42], [43]. The phase velocity of a beat wave can be 
modified in a natural way, by varying the frequency of the beating lasers (the use of two counterpropagating laser 
beams with a chirp has already been considered for the electron injection in a LWFA [44], but relying on a different 
mechanism than the one explored here). 

This Paper shows that, even at nonrelativistic radiation intensities (/ <C 10 24 W/cm 2 , in the case of IR lasers), 
ions can be accelerated to high energies directly by the ponderomotive force of an EM beat wave [42], [43] with 
variable phase velocity, driven by two counterpropagating lasers with variable frequency [45]. This scheme for 
acceleration in vacuum does not rely on complex laser-matter couplings and works in a single-particle regime, 
as in conventional accelerators, offering the possibility of direct, efficient control over the relevant beam features 
(e.g., mean energy, energy spread, and total accelerated charge) via regulation of the fundamental parameters of 
the two lasers (i.e., intensity, duration, and frequency variation). The technique can operate in a wide frequency 
range, including the IR and visible region of the EM spectrum, where extreme intensities are available, and it can 
be adapted to any source of ions: it can serve as an energy booster, accelerating ions from a beam with a given 
energy, or it can accelerate ions starting at rest from a tenuous gas or vapor that is turned into a plasma by the 
laser radiation. 

The Paper is organized as follows: in Section [II] the acceleration scheme is explored in detail with a one- 
dimensional theoretical model; in Section [Till the scaling laws for the relevant ion beam properties are presented 
and the requirements in terms of the laser features are discussed; in Section lTVl numerical results from self-consistent 
simulations in two dimensions are presented. Finally, in Section [V] the conclusions are stated. 

II. Theory of direct ion acceleration 

The acceleration technique proposed here is analyzed in the framework of a relativistic, one-dimensional, single- 
particle theory, in which the amplitude of both lasers is constant, referring to the configuration sketched in Fig. (Q~|). 
When considering the typical conditions required for the method to be effective, this model provides an accurate 
description of the fundamental features of the acceleration mechanism, as also confirmed by two-dimensional (2D) 
particle-in-cell (PIC) simulations (cf. Section llVT i. 

A. Field configuration and basic equations 

In the present theory, the lasers are described as two coherent EM waves propagating along the x direction, with 
vector potentials 

Aj = Aj{ cos (0) sin (£,-)] e y + sin (0) cos (&)] e^}, (1) 

where e y and e z are unit vectors, determines the polarization type (e.g., = ir/2 for linear polarization in the 
z direction, = 7r/4 for circular polarization), j = 1,2 labels the lasers, £i = x — ct (laser 1 propagates from 
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Fig. 1. (Color online) Sketch of the acceleration technique: a test ion is trapped and accelerated by the slow ponderomotive beat wave generated 
by two counterpropagating, variable-frequency lasers. Thin lines represent the vector potential of the two lasers (solid lines) and the sum of the 
vector potentials in the superposition region (dashed line); the thick, solid line represents the ponderomotive beat-wave potential. 



left to right), £2 = — x — ct (laser 2 propagates from right to left), Aj is the amplitude, and $ 3 are two 
arbitrary functions with first derivative 3^ > (henceforth, for the purpose of readability, the arguments £j 
will be omitted unless necessary). The local wavenumbers, kj (x,t) — ^f, and frequencies, uij (x,t) — ~^§f-, are 
k\ (x, t) = Ui (x, t) = c§[, fc 2 (x, t) = — <I> 2 , and w 2 (x, t) = c$ 2 - The superposition of the two EM waves can 
be described in terms of the two corresponding beat waves: (i) a fast, superluminal beat wave, having wavenumber 
K{x,t) = - $ 2 ), frequency Q(x,t) = §($i + $ 2 )> and phase velocity V^,(x,t) = § = c§f±|i, and (ii) 

a slow, subluminal beat wave, having wavenumber k(x,t) = h(G>i + 3» 2 )> frequency u>(x,t) = |($x — "3? 2 ) and 
phase velocity v ( j,(x,t) = ^ = c + ^, 2 (the following relations hold: kc — f2, Kc = u>, and v^V^ = c 2 ). The 
ratio between the frequencies of the two laser beams in x = at t = is assumed to be such that the phase 
velocity of the slow beat wave matches the initial velocity of the ion, c(3q, i.e., such that ^(0,0) = c/?o, leading 



to 



^2(0,0) _ l-Po 
wi(0,0) 1+Po ■ 

By resorting to the conservation of the transverse canonical momentum, p y = — - (Ai + A 2 ) • e y and p z 
I (Ai + A 2 ) • e z , the x component of the equation of motion for the ion is written, in cgs units, as 



dp x q 2 d 



dt 2 7 Mc 2 dx 



(Ax+A 2 ) 2 , (2) 



where q and M are the ion charge and mass, respectively, and 7 = y 1 + jji^s + q ^j^ 2 " 1 is the Lorentz factor. 
As shown in the Appendix, if the radiation intensity is not so high as to drive transverse oscillations with relativistic 
velocities (i.e., if Aj = j-^ <C 1), and as long as the characteristic time scales of the two beat waves can be 
separated, Eq. (f2]i can be averaged over the fast time scale to yield the ponderomotive equation of motion 

dp x iii 2 d . . 

— — = — COS ($1 — $2) 

dt 27 dx v ; 

= ^^fc(fi,t)Bill($i-*2), (3) 

7 

where dimensionless quantities (denoted with hatted symbols henceforth) have been adopted: k = with fc = 
k (0, 0), t = kgct, x = kox, and p x = Jj- [in these units, the initial frequencies of the lasers are cJi(0, 0) = 1 + /3o 
and W2(0,0) = 1 - (3 ]. In Eq. OJ, the averaged Lorentz factor is defined as 7 2 = 1 + p 2 x + A\/2 + A 2 ./2 + 
AiA 2 cos($i -$2). Equation Q can be obtained from the Hamiltonian Jf (x,p x , t) = [l +p 2 x + A\/2 + A\/2 + 
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A\a\ cos ($1 - $2)1 1 ^ 2 , using = — (the other Hamilton equation being % = %|- = %■); accordingly, 

the variation of the ion energy, -£L = ^f, is given by 

d 7 AxA 2 d . , 

— — = ^cos($i - <& 2 ) 

dt 2j dt V ' 

= -i^(&(5 ) ^)Biii($i-$a). (4) 

7 

Equation shows that the presence of a frequency variation is a necessary condition to effectively transfer energy 
from the beating EM waves to the ion over long times (cf., the Lawson-Woodward theorem [1], [46], [47]). If the 
frequency of both lasers is constant (i.e., if $' x and $ 2 are constants), Eq. dU reduces to 

dr = sin [ 2 ( 5 _ p t) + O ] , (5) 

dt 7 

where <f>o — $1 (0) — $2 (0) is the initial phase of the beat wave. In the reference frame where /3o = 0, Jf is time- 
indepedent, p 2 , + A\A 2 cos (2x + <po) is constantly equal to A\A 2 cos (<fio) (i.e., variations in p 2 are compensated 
by variations in p 2 + p 2 ), p x oscillates between ±P, with P = {A\A 2 [\ + cos^o)]} 1 / 2 , and 7 remains constantly 
equal to V = [1 + A\/2 + A 2 1 /2 + A\A 2 cos^o)] 1 / 2 ; in any other frame, .^f is time-depedent, p x oscillates between 
(PqT ± P)/(l — /3q), and 7 oscillates between (T ± 0qP) /(l — /?□)■ Thus, 7 cannot grow over long times. However, 
if at least one of the lasers has variable frequency, causing the phase velocity of the ponderomotive beat wave to 
vary, Jrf? is time-depedent in any reference frame, and continuous energy transfer from the beating EM waves to 
the particle is possible. For appropriate choices of $1 and $2, the beat wave can trap the ion and accelerate it to 
large energies. 

B. Resonant solutions 

A resonant solution of Eq. (0, X (t) , is defined by the resonance condition 



X (t) - t 



-x (i) - i 



(6) 



which imposes exact phase-locking between the ion and the ponderomotive beat wave. In order to find a particular 
resonant solution, one can choose the expression for $! (or $ 2 ) arbitrarily and then seek the expression of <J? 2 (or 
<!>!) for which Eq. (0 and condition (O are simultaneously satisfied. By making use of (O, Eqs. (0 and (0 can 
be written, as 

^ = ^ S in(^o)($i + $ 2 ) (7) 

and 

jj£ = ^Bm (AO (*!-*£), (8) 

where -f 2 ± = 1 + Aj/2 + A\j2 + A\A 2 cos ((j> ), 7 n = 7/70-L = (l - Pi) , r y = ■yo±r (with r being the proper 



time, such that ^ — 7), and p y — p x /lo± = y 7 2 — 1- Assuming, without loss of generality, that $2 is known, 
Eqs. and ([8]) can be combined to obtain 

d 2 f r i 

^=^[-ei-2«(r„)J, (9) 
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where (1q = Al ^ 2 sin (Sq) and £1 = X — i. By integrating Eq. Q one can determine the resonant ion trajectory 

To J 



MA 2 

i ox 

in the parametric form (jt (rn) , t (Ty)| and then use Eq. © to determine $1. For a generic shape of $2, this 
process usually requires a numerical approach, because of the intricate dependence of t on Ty through the integral 
t = Jq 1 "" 7n (t||) dr,(. In the special case in which the frequency of laser 2 is constant (as in [45]), i.e., <&' 2 — 1 — /3o, 
the right-hand side of Eq. (O does not depend on t and the calculation can be performed analytically. Integration of 



Eq. © yields |i (t,,) = §■ (/xr,, - 2) r,,, where /Lt = Mo/7oy, with 7 „ = V 1 ~ and c o = Toy (1 - A)) = yT$> 
hence, p y (T y ) and -y M (Ty) can be expressed as 

A d ^ l-^(l-/JTy) 2 
Pll ( r ll) = "7" = — ^ 7; s— . 

dry 2c (l-/xryJ 

di l + cjj(l- M Ty) 2 
dry 2co(l-/xry) 

and, consequently, X (ry ) and £ (Ty ) are 

X (Ty) = ^- [eg (/iTy ~ 2) M Ty - 2 bg (1 - M Ty )] , (12) 
t (Ty) = - J- [eg (^Ty - 2) /iTy + 2 log (1 - /iTy)] . (13) 



Finally, by noticing that 1 — jt/Ty = y 1 + -X" and t can be written as functions of £1 and replaced in Eq. ©, 

yielding 



$1 (li) = + 7^ H [l + 2^o (1 + #>) | 



(14) 



It is useful to analyze the behavior of the resonant solution, given by Eqs. ([T0]i-(fT3l. in the early stage of 
the acceleration, t <C H" 1 , and in the asymptotic limit, t 3> pT 1 . For i <C pT l , Eq. (fT3l gives, to second 
order in i, Ty (f) — -J^riit 2 and the explicit dependence of py, 7y, and Xy on i is written as py (t) w 

7onA) + + 2^o7^ 2 ^ 2 ' 7n (*) ~ 7o ii + A)A** + 1+ 2^°'| | /3 ° M 2 ^ and (*) ~ A)* + ^^M* 2 - In the nonrelativistic limit 
(A) — > 0, 7oy — > 1), the ion accelerates uniformly and its energy grows quadratically in time as ^^i 2 . For t 3> A* -1 , 
the behavior of the solution depends on whether laser 1 and the ion are copropagating (/i > 0) or counterpropagating 



(/i < 0): if fj, > 0, Eq. ( [T3l gives, to leading order in t, Ty (t) « i 



1 — exp I — 2co/i£ 



and, consequently, 



* (*) ~ * - § and 7|| (*) « Pn (*) « 2^7 ex P ( 2c oMt- 4); if M < 0, Eq. (13) yields t„ (t) « 2^/^ and, 
consequently, X (f) w — i and 7y (t) w p y (f) « ij co\p\i. The behavior is asymmetric because, if /i > 0, the 
frequency must be increased to maintain phase-locking, causing a continuous increase of the ponderomotive force, 
whereas, if [i < 0, the frequency must be decreased, causing a continuous decrease of the ponderomotive force. 



C. Trapping 

Effective ion acceleration can occur also when $i and $2 do not match the resonant solution (exact phase- 
locking). In fact, it is sufficient to choose $1 and $2 appropriately in order to guarantee that the ion trajectory stays 
close to the beat-wave trajectory until reaching the desired energy, i.e., that the ion is trapped by the beat wave. 
For given $1 and $2, the beat-wave trajectory x^, m (namely, the trajectory of the point of the beat wave having 
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constant phase </>o) is determined by solving the equation <I>i (x$ — t) — $2 (— x^ —t) = 4>o with respect to x<p 
[for resonant trajectories, X m = x^ (£)]. In order to determine the conditions for the occurrence of trapping, 
Eq. ^ is expressed in terms of the phase difference between particle and wave, ijj = 2 [x — x^, m], and using 
the proper time r, as 



d 2 tp 
d^ 



d_ 



U(1>,t), 



in which 



U(ip,T) = 2A 1 A 2 cos{ $1 



Cl,0o ( T ) + 77 



+ 2a <Po ( T ) ^. 



(15) 



(16) 



where £ Mo , £ 2 ,0 O ' and a 4> depend on r through i as £ 10o (r) = x 0o (i) -i, £ 2 ,0 O (r) = -i^, (f) and a 0o (r) = 
^3-^0 (*) = 7 (*) 57 7 (*) ^^0o (*) ■ According to Eq. (TTBT l. trapping is allowed only if the frequency variation 
is slow enough to guarantee that the effective potential U presents local minima. For trapped ions, \ip\ -C 2^^ 
and U can be approximated as 



U (ip, t) w 2iii 2 cos |fc [x^ (i) , i] ip + O } + 2o^ (I) V 
leading to the necessary condition for ion trapping 

|a0 o (t) I < a M (i) = A x A 2 k [x^ (t) , t\ , 



(17) 



(18) 



where a M (t), is the maximum value of the ponderomotive force. Thus, trapping is allowed only if the maximum 
ponderomotive force of the beat wave is greater than the inertial force associated with the beat-wave acceleration. 
The particular situation in which | a$ (i) | = a„ (i) corresponds to a resonant solution in the particular case 
<f>o = ±tt/2. When $1 and $2 correspond to a resonant solution and |</>g| ^ ir/2, regions where trapping is 
possible always exist, although the solution is stable for cos(0o) < and unstable for cos(</>o) > 0, with ip = 
corresponding to the bottom and top of a potential well, respectively. The relative width of the potential well 
in U (ip, 0), Aip, can be used to estimate the trapping efficiency r/a (defined as the fraction of trapped ions) as 
t] a - « — 1 — f: arcsin[a0 o (0) /a M ]. Examples of typical trapped and untrapped trajectories are shown in Fig. |2] 



D. Ion acceleration with linearly chirped lasers 

An important case of variable-frequency lasers, particularly relevant in experiments, is that of linearly chirped 
beams, in which $1 and <J> 2 take the form $1 U^j = O i + (1 + Po) €1 + "ill and $ 2 = 002 + (1 - Po) 6 + 
o~2%, where 0oj are constant phases and Oj are the chirp coefficients. The local frequencies and wavenumbers are 
given by Cj\ fx, t) = k\ (x, t) = 1 + po + 2a±^i and Cj 2 (x, t) = —k 2 (x, t) = 1 — Po + 2a 2 £ 2 . The wavenumber and 
the frequency of the slow beat wave are k (x, i) — \+a_x — a + i and u> (x, t) = Po+<r + x — a_i, where o_ = o\ —o 2 
and <t + = a 1 +a 2 . The phase of the ponderomotive beat wave is $1 — $2 = 0o + 2 (x — pot — a + xt) +a_ (x 2 + t 2 ), 



with 0o = 0oi — 002- Consequently, the beat-wave trajectory can be written as x<^ (t) as 

1 



/-\ cr + - 1 
^0o (*) = — * - ~ 



l-2(a + ~p a-)t+(a 2 



(19) 
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Fig. 2. (Color online) Space-time evolution of the dimensionless ponderomotive force — cos (<l>i — $2) (color scale, red and blue denoting 
accelerating and decelerating regions, respectively), and trajectories of a resonant ion (solid line), a non-resonant, trapped ion (green markers), 
and a non-resonant, untrapped ion (yellow markers), for a situation in which (3q = 0, $2 = 1> an & ^1 is given by expression 1141 . with 
A1A2 = 5 X 10 — 4 and </>o = 5tt/6. The inset represents the effective potential U at t = 0, gray marking the trapping region, and the 
corresponding initial positions of the three ions. 



If only one laser is chirped, i.e., if 02 = or u\ = 0, Eq. $1% reduces, respectively, to 



x { } ] (i) = 



* {2) d) = -i+ — 



t — 



(20) 
(21) 



l-yi-2(r a (l + /Jb)i 

Equations (|T9l)-(l2TI) suggest that the acceleration can be improved by chirping both lasers, with o\Oi < 0. When 
doing so, the specific values of o\ and 02 must be chosen with some care because the chirp in the two lasers 
affects the trapping efficiency in a different way. This is readily seen by writing the inertial term a<f, m in the 
case of ions trajectories close to xP (t) and x? (t): for ions traveling with approximately the beat-wave velocity, 



«0o (*) ~ 70o (*) I po (*) ^*0o © 

i^ , yielding 



where 7^ is the Lorentz factor associated with the beat-wave trajectory 



a 



-a n /l-2 ( T 1 (l-/3o)t 



/3 - 1 + 2^1-2(71(1- A)) « 



(2) 



<T2\/l-2tT 2 (l+/3 )t 



(22) 



(23) 



l + A)-2yi-2<r 2 (l+A)* 

According to Eqs. (l22b and ( f23l , when the ion and the chirped laser are copropagating (e.g., ax < and 0-2 = 0), 
a0 o (t) is monotonically decreasing in time, causing the trapping regions to widen; on the contrary, when the ion 
and the chirped laser are counterpropagating (e.g., o\ — and 02 > 0), a,p (t) is monotonically increasing in 
time, causing the trapping regions to narrow. When chirping both lasers {cr\ ^ and 02 7^ 0), the use of a stronger 
chirp in the laser copropagating with the ion is in general convenient: in such a situation, the general expression for 
ct<p (i) has a minimum for a given time (which is t = if <j\ = —02), such that the width of the trapping regions 
increases until m reaches its minimum, and starts decreasing immediately after. This is illustrated in Fig. [3] 
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showing the amplitude of the trapping region, Aip, as a function of time when only one laser is chirped (case A) and 
when both lasers are chirped (case B). In case A, Atp grows monotonically towards its maximum value, Atp = 2ir. 
In case B, Aijj reaches a maximum and then starts decreasing; when A^ vanishes the ion loses the beat wave and 
the acceleration stops. The corresponding evolution of the ion energy as a function of the acceleration distance is 
shown in Fig. [4] in case A, the ion acceleration continues indefinitely, but with decreasing accelerating gradient; in 
case B, the ion accelerates with approximately constant accelerating gradient until losing the synchronization with 
the beat wave. 




2 4 6 8 10 12 
Time [ps] 



Fig. 3. (Color online) Width of trapping region as a function of time when only laser one is chirped with u\ = —A\Ai (case A, blue/dark 
curve), and when both lasers are chirped with a\ = — O.65A1A2 and <T2 = 0.35j4ij42 (case B, red/light curve). Ion type: proton; I\ = I2 = 
10 20 W/cm 2 . 




0.2 0.4 
Distance [mm] 



Fig. 4. (Color online) Evolution of ion energy as a function of the acceleration distance for the same configurations, case A and case B, as 
in Fig. f5] Circles indicate equal time instants for both cases. 



The maximum energy gain, Ae M , occurs when the maximum chirp that allows for trapping is used, hence, accord- 
ing to Eq. <fT~9b , when \a_\ = A\Ai. For non-relativistic ions, the maximum energy is then Ae M as | fjjj^o) ~ 
i \ A1A2tj , leading to the scaling law 

Ae M [MeV] w 0.8 Z^Ar 3 ^ [10 20 W/cm 2 ] I 2 [10 20 W/cm 2 ] 

x A01 [Atm] A02 [/H AT 2 [ps] , (24) 
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where A\ is the ion mass number and Acy = 2tt /\kj ; (0, 0) are the initial laser wavelengths. Similarly, the 
acceleration distance, Ax k, |x0 o | ~ \A\A2i 2 , scales as 

Ax [ M m] « QZfArHl' 2 [lO 20 W/cm 2 ] 7 2 1/2 [lO 20 W/cm 2 ] 

x A01 [Mm] A 02 [H AT 2 [ps] . (25) 

By varying o\ and era, the energy gain Ae can be controlled, according to Ae = J Ae M . When increasing 

|er_|, a0 o increases and, consequently, the trapping efficiency decreases; hence, the chirp coefficients must be 
tuned in order to obtain the desired trade-off between energy gain and total accelerated charge. 

III. Beam properties and laser requirements 

By analyzing the range of validity of the present model, conclusions can be drawn for more general situations, 
in which spatial distributions of ions are considered and the effects due to the finite-size of the laser pulses are 
taken into account. 

In the presence of a spatial distribution of ions, the scheme remains effective if the space-charge field due to the 
ion density ri\ does not perturb the EM beat-wave structure. Assuming that no neutralizing electron background is 
present, this establishes the criterion 

rn [10 19 cm" 3 ] « Zr^ll' 2 [10 20 W/cm 2 ] 7 2 1/2 [lO 20 W/cm 2 ] . (26) 

When taking into account realistic laser beams, having finite duration Tj and spot size Woj, both the longitudinal 
and the transverse dependence of the laser envelope must be considered. As far as the longitudinal dependence is 
concerned, since long pulses (e.g., Tj > ps with visible/near-IR lasers) are required to obtain significant energy 
gains, the ponderomotive force associated with longitudinal intensity variations is negligible when compared with 
the ponderomotive force of the beat wave. As a result, A\ and A2 can be considered as slowly varying functions 
of t, and the only relevant effect is that intensities high enough for ion trapping exist only during a fraction of the 
overlap time for the two laser pulses; consequently, the duration of the acceleration process is limited to a fraction 
of the pulse duration, determined by the specific chirp law. As an effect of the transverse variations of the radiation 
intensities, the ions are pushed outward by the transverse ponderomotive force of the laser pulses, whose magnitude 
scales as Aj/Woj', hence, if the intensities of the two lasers are comparable, the ratio between the transverse forces 
and the ponderomotive force scales as 2 J^ ^ 1' w i m ^0 = 2-7r/fc (0, 0). In these conditions, the transverse cross 
section of the accelerating region is on the order of ttWq, where Wo is the overlapped spot size of the lasers, 
and the length of the acceleration process is always limited to distances on the order of z, = min (z r i, z^), where 
Zrj = irWQj/Aoj are the Rayleigh lengths of the laser pulses; therefore, the maximum total accelerated charge Q M 
(obtained when the laser duration is on the order of z r ) scales as 

Q M [pC] « SrfcZim [10 19 cm" 3 ] W% [/im] z r \pun] , (27) 

where 7] tr is the fraction of ions trapped by the beat wave (cf. Section 111-0 . 
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In light of these considerations, the scaling law (1241 can be conveniently expressed in terms of the laser-pulse 
energies as 

A£ M [MeV] « ,o 8 | ^goi^ 

ra o/JilWL. (28) 

A? z r i [/im] z r 2 [jjm] 

Hence, if Wjy/Aoj (and hence z r i) are fixed, the maximum energy gain depends only on the energy of the two 
laser pulses. For given pulse energies, wavelengths, and spot sizes, the combination of pulse lengths and intensities 
can be tuned in order to obtain the desired balance between the relative width of the final energy spectrum (which 
scales as 4^|f) an d tne accelerated charge (which scales as AiA 2 if the laser durations are longer than z r , 
being constant otherwise). Furthermore, the scalings indicate that, for given pulse energies, and keeping Woj/Xoj 
constant, higher energies are obtained when using lasers having shorter wavelengths. As an example, with 800 nm 
lasers (e.g., TkSapphire lasers) focused at 10 /im, an energy gain A£m = 1 MeV (with protons) requires 440 J in 
each laser pulse; with 10 /im lasers (e.g., CO2 lasers) and with the same Woj/Xoj, the same energy gain requires 
5.5 kJ in each laser pulse. 

Provided that the available laser energy is sufficient, the other fundamental requirement of the acceleration method 

is the frequency bandwidth. The total frequency excursions, Auj, necessary to accelerate an ion from the initial 

velocity [3q to the final velocity (3 are determined by 

= 2/3 + Am - Am 2 
2 + Alux + Alu 2 

When only one laser is chirped, Aojj are given by Acji = 2 (/? — /3o) / (1 — /3) and Aui 2 = 2 (/? — /?o) / (1 + /?)■ 
The minimum frequency excursion per laser is obtained when Acji = — Au>2 = /3 — 0q. 

The above estimates indicate that suitably chirped pulses, with energy level between 100 J and 1 kJ and with 
frequency excursions on the order of 5-10% of the central frequency (potentially available with current laser 
technology), should already allow for proof-of-principle experiments producing monoenergetic ion beams in the 
100 keV - 1 MeV energy range. To achieve relativistic energies, highly energetic light sources capable of delivering 
coherent radiation over a large frequency range are required. Alternatives can be pursued, for instance, with multi- 
stage acceleration schemes. 

IV. Numerical results 

The effectiveness of the method and the validity of the estimates in Sec. [HI] have been confirmed with one- 
dimensional (ID) and 2D, EM PIC simulations, using the Osiris 2.0 simulation framework [48], accounting self- 
consistently for space-charge and propagation effects. Here, a selection of results is presented from a 2D simulation, 
in which the source of ions to be accelerated is provided simply by a slab of hydrogen, with a thickness of 75 /im 
and proton density n\ = 5 x 10 16 cm~ 3 (in the simulation, field ionization of the neutral hydrogen atoms is calculated 
by using the Ammosov-Delone-Krainov rates [49]). Two approximately gaussian laser pulses are employed, with 
spot sizes Wqi = W02 = 10 /im and durations T\ = T 2 = 4.2 ps, focused on the mid plane of the slab. Laser 1, 
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propagating from left to right, has peak intensity I\ = 1.3 x 10 21 W/cm 2 and it is linearly chirped according to 
^1 (£i) — ^*oi£i + — £oi) 2 > where fcoi = 27r/Aoi, with Aoi = 800 nm, where £i indicates the distance from 

the pulse center, o\ = — 2 x 10 -5 fc^ is the chirp coefficient, and £oi = 6.25 x 10 3 fc^ 1 is the point of the pulse 
where k\ = fcoi, which needs to be specified when considering finite length pulses. Laser 2, propagating from right 
to left, has peak intensity I2 = 8.5 x 10 20 W/cm 2 , constant frequency, and wavelength A2 = 800 nm. The two laser 
pulses are launched from opposite x boundaries of the computational domain and hit the left face of the hydrogen 
slab at the same time. The computational domain is rectangular [L x = 300 /im (5 x 10 4 cells) and L y = 50 /im 
(100 cells)], with open-space boundary conditions in x and periodic boundary conditions in y. The left face of the 
hydrogen slab is located at x — 40 /mi. 




50 100 150 200 250 
x [/im] 



Fig. 5. (Color online) (a) Ion distribution in the x-p x phase space and (b) mean kinetic energy of the ions as a function of position, after 6 ps. 

The trapping process begins after 1.8 ps, with 4.2% of the protons initially within the focal region of the lasers 
being extracted from the bulk distribution and accelerated by the ponderomotive beat wave. As the acceleration 
continues, the trapped protons form a bunch that separates from the main distribution both in space and in 
momentum, as illustrated in Fig. After 6 ps, the proton beam has a mean energy of 12 MeV with a 7.5% 
energy spread (cf., Fig. [6] showing the proton spectrum), approximately 60 /im long and 20 /im wide, with a 
transverse momentum spread on the order of 10~ 2 Mc. Much lower energy spreads, down to 1% and below, can 
be obtained by employing longer pulses with the same energy and lower intensities (as discussed in Section [TTTb. 
The expansion of the remaining plasma bulk, visible in Fig. [5^, is driven by the space charge that forms because 
of the rapid expansion of the hot electron population produced by the lasers (in the range of radiation intensity 
considered here, the electron dynamics in the beat wave is dominated by stochastic heating [50], [51]). These results 
confirm the potential of the proposed technique for direct acceleration as a method for producing high-quality ion 
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Fig. 6. Distribution of proton kinetic energy after 6 ps. 

beams in a controlled way. Other simulation results, to be object of future work, suggest that, under appropriate 
conditions, higher plasma densities, generating significant space-charge fields, can even improve both trapping and 
acceleration. 

V. Conclusions 

The theory and the results presented here demonstrate the possibility of direct ion acceleration by the superimposed 
EM field of two counterpropagating, variable-frequency lasers with nonrelativistic radiation intensities. Calculations 
and simulations indicate that proof-of-principle experiments could already be performed using current ultraintense 
laser technology. The essential simplicity and robustness of the underlying physical mechanism, as well as its unique 
control capabilities, makes the present acceleration technique a promising candidate for the future development of 
laser-based ion accelerators, particularly when fine control over the beam features is required. 
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Appendix 

Validity of the ponderomotive approximation 
The squared vector potential, A 2 = (Ai + A2) 2 , is given by 



2 _ A\+A\ , A A , , sin 2 ((?)-cos 2 (c?) 



A = 



A Y A 2 cos($ s ) 



2 

[A\ + A 2 2 ) cos ($ s ) cos (#p) - {A\ - A 2 2 ) sin ($ s ) sin (3> F ) 



2A X A 2 cos ($ F ) 



(30) 



where <I>f(2;, t) = $1 + $2 and &s(%, t) = $1 — $2 indicate the phase of the fast and slow beat waves, respectively. 
In order to obtain the ponderomotive equation (01, the time average of the fast terms, containing $f, in Eq. d30l 
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must be negligible, such that one can assume A 2 = | (A 2 + A 2 ) + A1A2 cos($s) [which is rigorously valid only 
for circular polarization, i.e., when sin 2 (9) = cos 2 (9)]. This occurs if the two following conditions are met: (i) the 
variation of x and p x associated with the transverse quiver motion is negligible, which is verified if Aj = 4r-% <C 1 
(here, Aj = ^p-aj, where a = -f^ denotes the usual normalized vector potential, e and m being the electron 
charge and mass, and Z\ is the ion charge state); (ii) along the ion trajectory xj(£), the frequency of the slow beat 
wave is much lower than the frequency of the fast beat wave, i.e., lo\ = — | [^^-(ccj, t) + Vi^£-(xi, i)] is much lower 
than O; = — I [^ E (cci,£) + Wi^-(x;,t)], with Vi(t) being the ion velocity. The latter condition is automatically 
satisfied for trapped ions, because their velocity stays always close to the phase velocity of the slow beat wave, 
namely, Vi(t) v^[x\(t),t], leading to 

n i ~ c ($i + $' a )-^($' 1 -*' a )- ^ ; 

By recalling that v^, = c ^} + ^f , one concludes that ^ w for any t; therefore, Eq. <£3j provides an accurate 
description of the ion-trapping process, as well as of the long-term dynamics of trapped ions. 
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